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ANALYSIS AND EXPERIMENTS WITH A PULSED NEUTRON SOURCE FOR 
AN UNREFLECTED SOLUTION REACTOR UP TO $50 SHUTDOWN 
by Daniel Fieno, Thomas A. Fox, Robert A. Mueller, and C. Hubbard Ford 

Lewis Research Center 

SUMMARY 

Pulsed-source experiments have been performed using the NASA zero-power solu- 
tion reactor. This unreflected, homogeneous, thermal reactor provides a simple geom- 
etry for assessing the difficulties associated with obtaining shutdown reactivities by using 
the inhour method, the method of Simmons and King, and the various TT area-ratio Tt tech- 
niques. A large range of shutdown reactivity from $0 to $50 was studied. 

This study showed that multigroup- transport-theory- calculated values of the funda- 
mental prompt-mode decay constant agree with experimentally determined values to 
within 5 percent from delayed critical to $5 shutdown; from $5 to $50 the calculated 
values exceed experiments progressively from 5 to 15 percent. The inhour method and 
the fT area- ratio M techniques of Garelis and Russell and Gozani for obtaining the experi- 
mental reactivity gave results that agreed with each other over the range studied. How- 
ever, the methods of Simmons and King and Sjostrand considerably underestimated shut- 
down reactivity over the large range studied. Of all the methods of analyzing the experi- 
ments, the inhour method was the most accurate and easy to apply. 


INTRODUCTION 

There are two basic methods for experimentally determining the reactivity of a sub- 
critical multiplying medium that apply to the pulsed-source technique. One of the meth- 
ods determines the reactivity from the prompt- and delayed- neutron responses and gives 
rise to the so-called "area-ratio" techniques. The experiments which use the delayed- 
neutron response as part of the technique for experimentally determining the subcritical 
reactivity of a system are also called modified pulsed-source experiments. Three area- 
ratio techniques for obtaining subcritical reactivities are considered in this study: 

(1) the area-ratio technique of Sjostrand (ref. 1), (2) the extrapolated area-ratio tech- 



nique of Gozani (ref. 2), and (3) the area-ratio technique of Garelis and Russell (ref. 3). 

The other method for experimentally determining the reactivity of a subcritical 
multiplying medium makes use of the measured fundamental prompt-mode decay con- 
stant. Two techniques based on this method are considered in this study. In one tech- 
nique the subcritical reactivity is obtained by using the measured fundamental prompt- 
mode decay constant and a calculated parameter. This technique is referred to as the 
inhour method in this discussion and is described by Preskitt et al. (ref. 4). The other 
technique, that of Simmons and King (ref. 5), uses only the fundamental prompt- mode 
decay constants at the subcriticality in question and also at delayed critical. 

Modified pulsed-neutron-source techniques for experimentally determining the reac- 
tivity of a subcritical multiplying medium have been based on simplified reactor models. 
The area- ratio method of Sjostrand is based on monoenergetic diffusion theory applied 
to an unreflected reactor; the extrapolated area-ratio method of Gozani is based on the 
assumption that the spatial variation of the delayed modes is the same as that of the fun- 
damental prompt mode; and the technique of Garelis and Russell is also based on mono- 
energetic diffusion theory applied to an unreflected reactor. However, these experi- 
mental techniques are used for complex reactor systems which require correspondingly 
more complex reactor models. Attempts to improve the theoretical basis of the modi- 
fied pulsed-neutron-source technique have been made by Corngold (ref. 6), Becker and 
Quisenberry (ref. 7), Masters and Cady (ref. 8), Wallace et al. (ref. 9), and Preskitt 
et al. (ref. 4). 

The difficulties associated with the various modified pulsed-neutron-source tech- 
niques, as well as with the method of Simmons and King, indicate the need for data ob- 
tained for systems having a simple geometry. The NASA zero-power solution reactors 
are convenient for performing pulsed-source experiments and making corresponding 
calculations for systems ranging from delayed critical to highly subcritical. These re- 
actors use a solution of enriched (93. 2 percent uranium-235) uranyl fluoride salt 
(UOgFg) dissolved in water and contained in a cylindrical aluminum vessel. For a given 
concentration of uranyl fluoride salt in water, criticality is achieved by adjusting the 
height of the solution. No control rods are associated with these solution reactors. The 
geometry is thus ideal for performing pulsed-neutron experiments. The reactor used 
for the pulsed-neutron experiments reported herein is described by Fox et al. (ref. 10). 
In this report, various important aspects of the experimental procedure, the theoretical 
basis of the calculations, and the comparison of experimental results with calculations 
are discussed. 


DESCRIPTION OF EXPERIMENTS 


Pulsed-neutron-source experiments are straightforward to perform and give a max- 
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imum of information for the experimental effort involved. The equipment used for this 
study was high-quality, commercially available equipment with few problems of adapta- 
tion. Pulsed-source experiments require a great deal of care in their execution, since 
a number of variables must be considered. Considerable effort was required to deter- 
mine the operating variables such as pulse rate, pulse size, source and detector loca- 
tion, analyzer channel width, and operating time. It was also necessary to determine 
such inherent variables as reactor background, counting system resolution, and water 
evaporation from the reactor. 

In this section, the experimental arrangement, which varies somewhat in detail from 
the typical pulsed-source setup, is described. The experimental variables are then dis- 
cussed. Finally, a brief description of the experimental procedures and program is pre- 
sented. 


Experimental Arrangement 

Figure 1 is a block diagram of the general experimental arrangement used in this 
study. A pulsed-neutron source was located external to the solution reactor. Signals 
from a neutron detector inserted into the solution were sorted as a function of time by 
using a multichannel analyzer. The pulsed-source unit was used to trigger the analyzer 
sweep. This arrangement results in one analyzer sweep per pulse of the neutron gener- 
ator. The number of sweeps, as monitored by a scaler, therefore equaled the number of 
neutron pulses. 

Solution reactor . - Figure 2 shows an elevation view of the reactor tank with the 
pulsed-neutron source and the neutron detector in typical locations. The reactor is the 



Figure 1. - Diagram of pulsed-neutron-source experiments in NASA solution reactor. 
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Figure 2. - Typical pulsed-source geometry for unreflected NASA solution reactor. 
Drawn to scale. 


NASA Zero Power Reactor II (ZPR-II). The reactor tank consists of a 76.2-centimeter- 
diameter by 105- centime ter-tall cylindrical aluminum tank which is open at the top and 
surrounded by an annular region about 12. 5 centimeters thick. This annular region can 
be used to radially reflect the reactor with water, but no reflector was used in this study. 
The fuel is a water solution of UO 2 F 2 which can be varied over a wide range of concen- 
tration but, in this study, was kept fixed near a hydrogen to uranium-235 atom ratio 
(H/X) of 975. Reactor control is achieved by fuel manipulation. Therefore, a clean sys- 
tem and nearly ideal geometry are readily obtainable. A more complete description of 
the ZPR-II can be found in reference 10. 

Neutron pulse generator. - The neutron pulse generator is a small positive-ion 

” 9 A 

sealed-tube accelerator which produces 14. 3-MeV neutrons by means of the H (d, n)He 
reaction. The neutron output specifications include a maximum pulse size of the order 

7 

of 10 neutrons per pulse, a pulse width of 3. 5 microseconds at half amplitude, and a 
pulse rate variable from 1 to 11 pulses per second. Pulse size can be varied by chang- 
ing the acceleration voltage. Pulses are produced by feeding a voltage pulse to the step- 
up transformer which produces the ion accelerating potential. The active area of the 
target is about 5 square centimeters, which is small compared to the size of the reactor. 

Neutron detectors. - The detector probe which is inserted in the fuel solution is less 
than 1 centimeter in diameter and has a sensitive volume of 1 cubic centimeter. The 
small size approximates a point detector and introduces a negligible perturbation. The 
reactivity worth for this small probe has been measured as less than 7£ at the center of 
the core. The small size also makes it possible to find an optimum detector location 
which enhances the fundamental prompt mode in comparison with the higher-order modes. 
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The first small probe used in this study was filled with boron trifluoride (BFg) gas. 
It was found after a short period of use that tests at or near critical introduced flux 
fields of sufficient intensity to cause detector performance to deteriorate (the high- 
voltage plateau would deteriorate). This problem occurred with systems which were 
shutdown by less than $3 from delayed critical and got progressively worse as delayed 
critical was approached. 

O 

Substitution of an otherwise identical detector that is filled with helium gas (He ) 

was made. The data obtained were consistent and repeatable over the entire range of 

fuel solution heights including critical and near critical. For cases shutdown by $3 or 

3 

more, the data for the BFg and He detectors are very similar. Data from both probes 
are presented in this report. 

Multichannel analyzer . - The signals from the detector were fed through a preampli 
fier and amplifier to a 400- channel analyzer which was operated in the multiscaler mode 
In this mode, the channel width is variable. For these experiments, the channel width 
was varied between 50 and 400 microseconds with a fixed 12. 5-microsecond gap between 
channels to read the data into the memory. System resolution time was measured to be 
1. 12 microseconds, which made it possible to use a moderately large pulse of neutrons 
without encountering large deadtime counting losses. 


Experimental Variables 

The variables associated with this pulsed-neutron-source study are divided into two 
groups for discussion. Variables which can be readily controlled and changed are 
covered first. Then variables which are inherent to the physical setup are determined 
or accounted for. 

Experimental variables which can be controlled include neutron pulse size, pulse 
rate, number of neutron pulses, neutron source location, neutron detector type and de- 
tector location, analyzer channel width, and reactor solution height. Among the more 
important factors which can vary but which are not controlled by the experimenter are 
various backgrounds, room scatter, water evaporation correction, and neutron source 
changes with use and age. In addition, such factors as counting system resolution and 
analyzer channel gap must be determined. 

Neutron pulse generator variables . - The neutron generator has been described 
briefly in the previous section. The energy of the neutrons generated, the width of the 
pulse, the time delay between the electrical trigger pulse and the neutron pulse, and the 
physical size (target area) of the source are regarded as fixed. The pulse size, pulse 
rate, and total number of pulses are considered in this section. 

No attempt was made to measure the number of neutrons generated per pulse or to 
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determine quantitatively how much variation in pulse size occurred. The log N ampli- 
fier of the reactor control instrumentation was used to monitor the neutron level through- 
out the entire experiment. This instrument channel detected and displayed any average 
change in neutron level over periods of time that are very long compared to the time of a 
single pulse. The pulse size can be changed by adjusting the accelerator high-voltage 
pulse size from the control unit. This voltage setting was selected so that the deadtime 
count loss in the analyzer channel with the most counts was less than 10 percent. In ac- 
tual practice, only a few percent counting loss occurred in the highest count channels. 

The pulse rate is variable over a range from approximately 1 to 11 pulses per sec- 
ond. The effect of pulse rate on experimental results over a range of reactivities was 
determined. It was found that pulse rate had little effect on results for the range of 3 to 
11 pulses per second. From a practical consideration, the highest pulse rate possible 
should be used to minimize operating time. It was found, however, that 11 pulses per 
second apparently exceeded the duty cycle of the generator and caused the target to over- 
heat in a short period of time. It was necessary to reduce the pulse rate to 5 to 6 pulses 
per second in order to give maximum stability to the source output. 

The total number of pulses (the operating time) was determined by the requirements 
for good statistical data. This varies with the configuration involved. For configura- 
tions near delayed critical, the prompt -neutron response data determined the amount of 
pulsing required. For the highly shutdown cases, the delayed-neutron response was the 
limiting factor. The near critical cases do not present any difficulty since it was easy to 
get enough counts of both prompt and delayed neutrons in a short operating time. For the 
highly shutdown cases, it was necessary to pulse for as long as 90 minutes. In general, 
it was desirable to get about 1000 counts from delayed neutrons (no less than 2 counts in 
any channel) in order to have sufficiently good statistics for reactivity determination 
using any of the area- ratio methods. Adequate counts were obtained in the prompt 
region when this criterion was met in the delayed region for the shutdown cases. 

Analyzer channel width . - For the purposes of this study, the multichannel analyzer 
was always used in a multiscaler mode. The entire 400 channels were employed in all 
cases. Only the channel width, which was kept as wide as practicable in order to count 
as many neutrons as possible per pulse, was varied during this study. The channel 
width was adjusted in general so that at least 25 to 30 channels were available from 
which to determine the prompt- neutron response. 

Source and detector locations . - Preliminary tests were conducted with both the neu- 
tron source and a detector external to the reactor core. Although such tests provided 
useful information regarding the experimental arrangement, no satisfactory pulsed data 
were obtained. The data were found to be extremely sensitive to the source and detector 
locations relative to one another. The data were also affected by scattering of neutrons 
from the walls and structure in the room. It was necessary to insert the detector into 
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the solution reactor core to alleviate the room scattering problem. 

A limited study of the effect of detector position on pulsed -neutron results was car- 
ried out. The maximum effect on the decay constant as a result of detector position was 
of the order of 2 percent total variation. Since the reproducibility of the data is of a 
similar magnitude, a quantitative result is impractical. There appeared to be a consist- 
ent reduction in the decay constant as the detector was moved from the center to the edge 
of the core in a radial plane, or from midplane to top of the core along an axial line. 

This effect did not seem to be a function of the magnitude of reactivity shutdown. 

The effect of neutron detector location on the reactivity, as determined by the var- 
ious area-ratio methods, had too many inconsistent data points to ascertain a definite 
trend. As a consequence, the detector was located at the position which theory indicates 
should give best results. 

With the detector located in the solution, the effect of changing source location ap- 
pears comparable to changing the pulse size. Different source locations result in a dif- 
ferent number of neutrons being injected into the system, but the modal characteristics 
are essentially a function of system geometry rather than the location of the source. 

For these studies the solution had an H/X of 975, and the system had a length to 
diameter ratio of nearly 1/3. For such conditions the source was most effective when 
located on the axis of the cylinder near the reactor tank bottom. The detector location 
used was at a radius of about 16. 6 centimeters (0. 44 of the radius) and at a height about 
two-thirds of the core fuel height. The higher modes are minimal with respect to the 
fundamental mode at this location. 


Experimental Program and Procedures 

Using the experimental arrangement and experimental variable criteria described 
in the preceding sections, a series of experiments was made on reactor configurations 
ranging from the delayed critical condition to highly shutdown systems using a BFg de- 
tector. These experiments covered several days and were repeated in their entirety 

3 

several weeks later when a He detector was procured. 

When a solution reactor system such as this is operated with its top surface open to 
the atmosphere, there is a continual evaporation of water from this exposed surface. 

This evaporation results in a gradual increase in solution concentration which must be 
continually accounted for and periodically corrected for by the addition of water. This 
reactor change was determined by running a clean delayed critical case at the beginning 
of each day's experiments. Corrections to adjust the data to an H/X of 975 were made 
accordingly. This reference check also assured that the solution concentration was being 
kept close to an H/X of 975 to keep the corrections as small as possible. 
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Pulsed experiments a t delayed critical. - It was difficult to conduct pulsed-source 
experiments at delayed critical because of the constantly increasing neutron level with 
pulsing. It was desirable to achieve criticality at as low a neutron level as practicable 
for the initial condition. The prompt component was then more readily distinguishable 
from the steady- state neutron level when the system was pulsed. More important, if the 
initial neutron level is too high, the operating range of the neutron counting equipment 
may be exceeded and the deadtime counting losses become excessive before the experi- 
ment is completed. As an alternative, the technique suggested in reference 11 was tried. 
This involves a number of short runs of 30 seconds each in which the neutron level was 
reduced to some low level after each run. This method did not prove as satisfactory for 
our study as a single 5-minute run at about 5 to 6 pulses per second where pulse size was 
kept smaller than for the subcritical runs. 

Subcritical pulsed- neutron experiments. - The subcritical configurations were 
achieved by simply lowering the fuel solution level below the critical height. This re- 
sulted in an unperturbed core with a similar clean geometry as in the delayed critical 
cases. This arrangement permits a study of parameters as a function of fuel height (AH 
below critical) at as many different heights as desired. 

When a subcritical configuration is pulsed repetitively, an equilibrium level of de- 
layed neutrons is built up. This neutron level is a function of the pulse rate, the pulse 
size, and the subcritical state of the system. As soon as the pulsing ceases, this equi- 
librium level slowly decays to the reactor background level. 

In general, more configurations were run at conditions near delayed critical than at 
the highly shutdown state. Near delayed critical (delayed critical to-$5), the pulsing 
time required is short, and the backgrounds can be determined with greater ease and 
certainty. In order to obtain sufficient data for the highly shutdown cases, the running 
time per case becomes very long. 

When obtaining data to determine reactivity by any of the area-ratio methods, it is 
essential that both the prompt- and delayed-neutron ’’areas'' be obtained. The funda- 
mental prompt-mode decay constant and prompt area proved easier to measure than the 
delayed-neutron area. This is due, in part, to the poorer counting statistics for the 
delayed-neutron response. In addition, the decaying delayed-neutron contribution at the 
end of pulsing must also be included in the delayed data. The experimental arrangement 
employed in these tests made such a measurement of the postpulsing response difficult. 

It was determined by means of a series of experiments that prepulsing the system until 
its equilibrium delayed-neutron level was obtained (>3 min) provided a method that gave 
similar results without the need to measure the postpulse decay area. In addition, pre- 
pulsing was found desirable, if not necessary, to stabilize the pulse size. When the 
pulser is started after not pulsing for some time, it frequently takes time for the ion 
source to come into equilibrium . Until equilibrium occurs, the neutron pulse amplitude 
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can vary by orders of magnitude. Therefore, all subcritical data reported were obtained 
by using the prepulse technique. For the critical cases prepulsing was not feasible. 

Reactor bac kground . - For each subcritical configuration there is a neutron back- 
ground that is always present. The nearer to delayed critical, the higher this back- 
ground level will be. Background was measured at a fixed time after pulsing ceased. It 
has been determined that most of the cases studied were not particularly sensitive to the 
reactor background. The reactivity measurements did show some effect from an arbi- 
trary change in background introduced into the data reduction scheme. However, the 
fundamental prompt-mode decay constant did not change, even when reactor background 
was arbitrarily increased by a factor of 2 or more. 

For critical configurations a neutron level is achieved rather than a background. 
Because of the nature of critical systems, this can vary widely in magnitude depending 
on the steady-state power level present. Pulsed critical cases start with as low a level 
as practicable and end with a higher level, both of which are constant. This area, which 
corresponds to the delayed- neutron area in a subcritical case, is subtracted from the 
total data before the fundamental prompt mode decay constant is determined. 


ANALYSIS 

Analysis of Experimental Data 

For a given subcritical solution height, the reactor was pulsed repetitively until ade- 
quate counting statistics were obtained. Figure 3 is a schematic representation of the 
data obtained from an idealized pulsed-source experiment. The background-corrected 
logarithm of the detector response is plotted as a function of the analyzer channel num- 
ber. Shown is the equilibrium delayed-neutron contribution persisting after the decay of 
the prompt-neutron portion of the curve. The slope of the prompt portion of the curve 
yields the fundamental mode decay constant s y m t*ols are defined in appen- 

dix A. ) The areas corresponding to the time integrals of the prompt, as well as the de- 
layed, neutrons are also shown in figure 3. 

Appendix B indicates how the experimental data obtained for a subcritical reactor 
are related to reactivities. The reactivity in dollars for the area-ratio method of 
Sjostrand is obtained by taking the ratio of the prompt-neutron area to the delayed- 
neutron area (refs. 1 and 12); that is, 
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Figure 3. - Ideal detector response to pulsed neutron source. 
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The quantity Zq (F) is the spatially dependent kinetic distortion factor which is obtained 
by calculation using equation (B72')- The reactivity in dollars for the area-ratio method 
of Gozani is obtained by taking the ratio of the extrapolated fundamental prompt-mode 
area to the delayed-neutron area; that is, 
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The quantity Z^(F) is the spatially dependent kinetic distortion factor which is obtained 
by calculation using equation (B53). In the Garelis-Russell method, the prompt-neutron 
response is weighted by the factor exp(xt), and then integrated over time. The quantity 
x is adjusted so that this integral equals the prompt- plus delayed-neutron area; that is, 



N (F, t)exp(xt)dt = 

ir 
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Then the reactivity in dollars is given by 


P? R ($) 
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where the kinetic distortion factor Zq (?) is equal to Zq (r ). The quantity x in the 
Garelis-Russell method thus corresponds to the parameter /3q/Aq, the ratio of the ef- 
fective delayed-neutron fraction to the generation time for the fundamental prompt mode. 

The methods just described for obtaining reactivity are the so-called area-ratio 
methods. These area-ratio methods are also described as being modified pulsed- source 
experiments. Reactivity may also be determined from the inhour method by using the 
measured fundamental prompt-mode decay constant and a calculated value for the param- 
eter /3 q/Aq. This reactivity is given by 


Pg N ($) 



(B75) 


Preskitt et al. (ref. 4) discuss the merits of using the inhour method for determining the 
"experimental" reactivity. 

The experimental reactivity may also be determined by the method of Simmons and 
King (ref. 5). In this method, the fundamental prompt- mode decay constant must be 
measured at delayed critical and at the subcritical configuration in question. Then the 
reactivity is given by 


P® K ($) = - 



(B79) 


where 



is the fundamental prompt-mode decay constant at delayed critical. 


As indicated, appendix B shows how all of the various equations are derived, while 
appendix A lists and defines the symbols used. In particular, R is the repetition rate 
for the pulsed source. Also the kinetic distortion factor for the three area-ratio tech- 
niques is a quantitative measure of the spatial and energy differences between the 
delayed- and prompt- neutron modes. 

The reduction of the experimental data into the parameters of interest was carried 
out by using a digital computer. The computer program GRIPE II developed by Kaufman 
(ref. 13) was adapted for use on the IBM-7094. This program calculates the fundamental 
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prompt-mode decay constant, as well as the reactivity, by the area-ratio methods devel- 
oped by Gozani, Garelis and Russell, and Sjostrand as discussed previously. An error 
analysis of each parameter is also carried out by this code, as well as a quality check of 
the raw data. 

After the experimental data have been ordered according to time, they are resolu- 
tion corrected for counting losses. The reactor background was introduced as a constant 
for all channels. The experimenter makes a conservative estimate of the number of data 
channels over which delayed background appears to be constant. The code makes a sta- 
tistical search to find the best value of the delayed background. The average value of 
total reactor plus delayed-neutron background is then subtracted from each data point, 
thus giving the prompt-neutron data. These prompt-neutron data are next examined by 
the code to find the region of best fit to a negative exponential function. This exponential 
is assumed to be the fundamental prompt mode. The best fit is determined by a least- 
squares analysis of the experimental data in which the data are weighted in importance 
approximately according to the number of counts per data channel. The channels which 
determine the fundamental prompt-mode region of the decay curve are determined by an 
iterative calculation which successively rejects channels from each end of the exponen- 
tial portion of the curve until the best fit is obtained. The initial iteration is chosen to 
be conservatively large to assure that the best fit is clearly delineated. A chi-square 
goodness -of -fit test is then applied to assure that the data are described by an exponential. 

The reactivity values are determined for all iterations by the code. The reactivity 
values reported are for the iteration with the best decay curve fit. The error analysis is 
based on counting statistics only and is not indicative of the absolute accuracy of the ex- 
perimental data. 


Absolute Criticality Calculations 

The theoretical basis for the pulsed-source calculation is described fully in appen- 
dix B. The S n method as described in reference 14 was used to obtain solutions to 
equation (B24) for the fundamental prompt-mode flux and to equation (B13) for the static 
adjoint flux. One-dimensional slab calculations along the axis of these cylindrical reac- 
tor systems were performed by using this method. The aluminum bottom plate was ex- 
plicitly considered in the transport calculations. The approximation with elastic 
neutron scattering treated in the P-^ approximation was found to be adequate in treating 
the neutron leakage from the ends of these subcritical solution reactors. Radial leakage 
was calculated from the radial buckling of the system which incorporated transport- 
theory- calculated extrapolation lengths. 

These calculations used eight neutron energy groups, seven fast groups and one ther- 
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mal group which included an up-scattering transfer component. The fast group cross 

sections were obtained by using the GAM-II code (ref. 15), and the thermal group cross 

sections were obtained by using the GATHER- II code (ref. 16). Table I gives the energy 

boundaries for each neutron energy group. Calculations with this group split have been 

found to accurately predict the criticality of NASA solution reactors over a large range of 
235 

hydrogen to U atom ratios. 

Delayed-neutron data used for these calculations were obtained from Keepin (ref. 17, 
p. 90). Data for are listed in table II. For the fraction of delayed neutrons 

produced by fission /3 was taken to be 0. 0065. The average energy of the delayed neu- 
trons varies from about 0.25 to about 1 MeV (ref. 17, p. 95). Thus, all of the spectra 
for each delayed-neutron group were taken to be in energy group 3. 

The solution to equation (B24) gives the fundamental prompt-mode decay constant 
ajj. For a given subcritical system, equation (B24) is solved iteratively until an is 
found such that the reactor multiplication is 1. 0//3 . Then equation (B13) is solved for 

r 


TABLE I. - NEUTRON ENERGY GROUP SPLIT 


Group 

Energy 

Lethargy 

1 

14.92 MeV to 2.231 MeV j 

“0.4 to 1.5 

2 

2. 231 MeV to 0. 821 MeV 

1.5 to 2. 5 

3 

0. 821 MeV to 0.224 MeV j 

2.5 to 3.8 

4 

0.224 MeV to 9.12 keV 

3.8 to 7.0 

5 

9. 12 keV to 454 eV 

7.0 to 10.0 

6 

454 eV to 8. 32 eV 

10.0 to 14.0 

7 

8.32 eV to 0.414 eV 

14.0 to 17.0 

8 

0.414 eV to 0 eV 

17. 0 to 


TABLE II. - DELAYED-NEUTRON DATA 
FOR URANIUM-235 


[From ref. 17, p. 90; /3 = 0.0065.] 


Delayed 

group 

Decay 

constant, 

K 

sec* 1 

Relative 

abundance, 

Delayed- 

neutron 

fraction, 

% 

1 

0.0124 

0.033 

0.000 214 5 

2 

.0305 

.219 

. 001 423 5 

3 

. Ill 

. 196 

. 001 274 0 

4 

.301 

.395 

. 002 567 5 

5 

1. 14 

. 115 

. 000 747 5 

6 

3.01 

.042 

. 000 273 0 
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the static adjoint flux for the same subcritical system. This immediately gives for 

this system. The fluxes obtained by solving these two equations are then used to deter- 
mine the generation time Aq (by using eq. (B15a)) and the effective delayed-neutron 
fraction /3 q (by using eq. (B19)). These two calculations are sufficient to give the reac- 
tivity in dollars from either of the two equations 

K -I 

pC($) =- eff _ (B20a) 

K aff*0 

or 

a P 

Pg($)==^+ 1 (B20b) 

^0 


These two calculations yield the same result for the static reactivity and can thus be used 
to check the internal consistency of the computed parameters. 

In addition to the calculations described, a solution to equation (B12) for the static 
flux is needed if kinetic distortion is present in the modified pulsed- source experiments. 
This spatially dependent factor is given in appendix B for the reactivity methods of 
Gozani, Garelis and Russell, and Sjostrand. For these clean, homogeneous, unreflected 
solution reactors, the assumption is made that kinetic distortion is not present. 


RESULTS 

Calculated Parameters, 

Using the methods discussed, calculations and experiments have been performed for 

a range of fuel solution heights for the NASA Zero Power Solution Reactor (ZPR-II) at a 
235 

hydrogen to U atom ratio of 975. The calculated variation with fuel solution height of 
the reactor parameters required is presented in figure 4. Calculated values of the effec- 
tive multiplication factor are shown in figure 4(a). The varies from about 

0. 64 to 1. 00 for a fuel solution height varying from about 13 to 26. 8 centimeters. The 
calculations gave a fuel solution height of 26. 80 centimeters for the reactor at delayed 
critical. Figure 4(b) shows the fundamental prompt-mode generation time Aq as a 
function of fuel solution height. Over the range of fuel solution heights studied, the gen- 
eration time varies from about 92 to 64 microseconds. Figure 4(c) shows the variation 
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of the fundamental prompt-mode effective delayed- neutron fraction j3q with solution 
height. The effective delayed-neutron fraction /3 Q varies from about 0. 0082 to 0. 00736 
over the range of fuel solution heights studied. The value of /3 was taken as 0. 0065 for 
U . Figure 4(d) shows the variation of the parameter £q/Aq as a function of fuel so- 
lution height. This parameter varies from about 89 to 115 per second over the range of 
fuel solution heights studied. 


Comparison of Experiment and Calculation 

The experimental results obtained from the pulsed-source study are given as two 
parameters: (1) the fundamental prompt-mode decay constant and (2) shutdown reactivity 
values by several techniques. All experimental results are compared with calculated 
values. The kinetic distortion factor has been taken as 1 for the three area-ratio meth- 
ods for all subcritical configurations. 

Fundamental prompt- mode decay constant . - Figure 5 shows the calculated and ex- 
perimental values of the fundamental prompt-mode decay constant in seconds - '*' as a 
function of difference of solution height from delayed critical AH in centimeters. The 
parameter AH has been used rather than reactor height for convenience. The experi- 
mental critical height at an H/X of 975 has been measured in ZPR-II as 26. 94 centi- 
meters as compared to the calculated 26. 80 centimeters. 

Small variations in reactor temperature and solution concentration from one experi- 
mental run to the next are unavoidable with the ZPR-H reactor arrangement. In reducing 
the raw experimental data, the solution heights were first corrected to room tempera- 
ture, 20° C, by using the experimentally determined temperature coefficient reported in 
reference 10. The AH values for a given experiment were then normalized to an equiv- 
alent AH at H/X of 975 by the relation described in appendix B of reference 18. Both 
corrections are small since the temperature and concentration were not permitted to vary 
over a large range. These corrections, though small, reduce significantly the scatter 
present in the uncorrected experimental data. 

The experimental data points for in figure 5 follow a smooth curve and are in 
good agreement with calculations. The value of at delayed critical is -115 seconds - '*' 
for calculation as compared to an average experimental value of -116 seconds - '*'. 

Several data points were repeated a number of times during the course of this study 
for various reasons. These repeated points give an indication of the reproducibility of 
the data. An inspection of the data given in table HI shows that three data points were re- 
peated exactly. However, a number of points were close enough in AH that adjustments 
could be made to bring the points to the same AH for comparison. In the range of AH 
from about 1 to 10 centimeters, the standard deviation in is less than 1 percent of 
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Figure 5. - Fundamental prompt-mode decay constant as 
function of fuel solution height from criticality at 
hydrogen to uranium-235 atom ratio H/X of 975. 


the value of ajy Beyond 10 centimeters, the standard deviation increases slightly but 
does not exceed 1. 5 percent at the largest AH studied. The error limits shown in ta- 
ble III are the standard deviations caused by counting statistics and represent only part of 
the total problem of reproducibility of the data. For points less than 1- centimeter shut- 
down, except for the delayed critical case, no deviation determination is possible. In 
this range, the percentage of uncertainty in AH can become a significant factor in the 
reproducibility of a data point. 

A smooth curve drawn through the (Kq data of figure 5 for AH values less than 
1 centimeter, but excluding the delayed critical data, was extrapolated to a value of about 
. -115 seconds” 1 for at delayed critical. This is in good agreement with the measured 
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TABLE m. - EXPERIMENTAL AND CALCULATED VALUES OF 


FUNDAMENTAL PROMPT-MODE DECAY CONSTANT 
[Unreflected 76. 2-cm-diam core with fuel solution of H/X = 975. ] 


Fuel solu- 
tion height 
from 

criticality, a 
-AH, 
cm 

Fundamental prompt-mode 
decay constant 

Fuel solu- 
tion height 
from 

criticality, a 
-AH, 
cm 

Fundamental prompt- mode 
decay constant 

Experimental, b 
sec'^ 

Calculated, 

_ a P 

5 

sec 

Experimental, b 
sec” 1 

Calculated, 

a 0’ 

sec” 1 

c 0 

115±0. 6 

116 

6.61 

1904±5 

2080 

c 0 

117±0. 6 

116 

6.65 

1937±6 

2098 

C . 067 

126±0. 6 

130 

7.62 

2297±7 

2500 

c . 115 

136±0. 5 

141 

8.18 

2496±12 

2760 

c . 251 

166±0. 5 

172 

8.22 

2524±9 

2780 

c . 508 

218±0. 6 

230 

8.33 

2586±8 

2830 

1.00 

330±1 

345 

9.16 

2918±10 

3230 

c l. 03 

330±1 

350 

C 9. 17 

2938±8 

3240 

c l. 04 

335±1 

355 

9.20 

2950±11 

3250 

2. 04 

568±2 

600 

9.37 

3043±10 

3340 

c 2. 07 

572±1 

610 

10.18 

3417±29 

3800 

c 2. 11 

5 92 ±2 

617 

10.81 

3797±15 

4180 

c 2. 11 

591±2 

617 

10.98 

3828±15 

4290 

3.05 

815±2 

865 

10.98 

3820±10 

4290 

4.06 

1096±4 

1170 

11.77 

4252±15 

4830 


1092±3 

1170 

c 12. 20 

4607±19 

5120 

4.06 

1091±4 

1180 

12.22 

4604±16 

5140 

C 4. 08 

1404±6 

1520 

12.58 

4782±10 

5400 

6.10 

1732±5 

1880 

12.74 

5008±12 

5525 

6.18 

1749±6 

1910 

12.88 

4977±17 

5640 

C 6. 60 

1914±7 

2075 





Experimental delayed critical height, 26. 94 cm; calculated delayed critical height, 


26. 80 cm. 

^Experimental errors based on counting statistics only. 
°B.e^ detector data, all other data with BF^ detector. 



values of at delayed critical which can be found in table HI. 

The agreement between experimental and calculated ajj is best at delayed critical 
and poorest for the highly shutdown cases. The calculated shutdown values are always 
larger in absolute magnitude than the experimental values of for this fuel concentra- 
tion. The percentage of difference between calculated and experimental values is less 
than 10 percent for cases as highly shutdown as $30, or an of -3000. The difference 
is less than 15 percent out to $50 shutdown, or an of -5000. 

The fundamental prompt-mode decay constant is a very sensitive parameter espe- 
cially near delayed critical. Some idea of the sensitive nature of this parameter can be 
shown by comparing the difference in calculated and experimental values to the corre- 
sponding change in the multiplication constant K g ^. For example, near critical a 5 per- 
cent change in the decay constant represents about 0. 04 percent change in K g ^; at about 
1 centimeter shutdown a 4 percent change in the decay constant is about 0. 1 percent 
change in K and at about 12 centimeters shutdown a 13 percent change in the decay 
constant is about a 4. 8 percent change in K g ^. The numbers given at 1 and 12 centi- 
meters shutdown are the differences between experimental and calculated decay constants. 

Determination of s hutdown reactivity. - As discussed in the INTRODUCTION, two 
general methods of determining reactivity values from the experimental data have been 
used. One method uses the values of to determine the reactivity, and the other uses 
an area-ratio technique of one form or another. Altogether, five different determina- 
tions of shutdown reactivity worth are presented for each case measured. 

Figure 6(a) is a plot of reactivity as a function of the shutdown fuel height increment 
AH in which calculations with the inhour method and the method of Simmons and King are 
compared. The solid curve shows the calculated values while a dashed line is used to 
connect the symbols showing the experimental data. The best agreement with calcula- 
tions is given by the inhour method p^($), which uses the experimental and calcu- 
lated values of /3 q/Aq given in figure 4(d) to determine shutdown reactivity. The 
Simmons-King method p?^($) gives results in poorer agreement with calculation. This 
technique assumes that /?q/Aq remains constant for all shutdown cases. 

Reactivity data as derived from the various area-ratio methods are compared with 

calculated and inhour method data in figure 6(b). The data given by the Gozani method 

p o ($) are in the best agreement with the results given by the inhour method and calcula- 
s 

tions, while the Garelis and Russell pj 1 ($) data approach this agreement. The data 

SO fa 

from the Sjostrand method p„ ($) are in poor agreement with calculation but are in fair 
agreement with the Simmons and King data. The reactivity data from all cases consid- 
ered in this study are given in table IV. More experiments were carried out than are 
presented in figures 5 and 6. This larger number of experiments was considered neces- 
sary in order to establish a degree of confidence in the methods. 

In considering the reactivity values derived by using the area-ratio methods, it is 
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(b) Comparison of calculation with area-ratio methods. 


Figure 6. - Shutdown reactivity as function of solution height from criticality 
at hydrogen to uranium -235 atom ratio H/X of 975. 



TABLE IV. - EXPERIMENTAL AND CALCULATED VALUES OF SHUTDOWN REACTIVITY 
[Unreflected 76. 2-cm-diam core with fuel solution of H/X = 975.] 


Fuel solution height 



Shutdown reactivity in dollars 


from criticality, a 
-AH, 

Calculated, 

p 



Experimental 



cm 

-Pb 

Inhour 

method, 

Simmons -King 
method, 

Area-ratio methods 1 

D 








-P m 

H s 

-pf 

Sjostrand, 

-p so 

Garelis-Russell, 

-p® 

Gozani, 

_G0 

" p s 

c 0. 067 

0.13 

0.10 

0.09 

0. 14±0. 001 

0. 14±0. 001 

0. 14±0. 001 

c . 115 

.23 

.18 

.17 

. 21±0. 001 

. 21±0. 001 

.21±0. 001 

c . 251 

.50 

.45 

.43 

. 44±0. 001 

. 46±0. 001 

. 46±0. 001 

c . 508 

1.01 

.90 

.88 

. 95 ±0. 004 

1. 00±0. 004 

1. 00±0. 004 

1.00 

2.03 

1.89 

1.84 

1. 8l±0. 01 

1. 90±0. 01 

1. 90±0. 01 

2.04 

4.31 

4. 03 

3.90 

3. 79±0. 05 

4. 03 ±0. 05 

4. 04±0. 05 

c 2. 11 

4. 46 

4.24 

4.10 

3. 96±0. 05 

4. 22 ±0. 05 

4. 23 ±0. 05 

2.11 

4.46 

4.23 

4. 09 

3. 83 ±0. 04 

4. 07 ±0. 04 

4. 08±0. 04 

3.05 

6.72 

6.28 

6.03 

5. 77±0. 12 

6. 20±0. 12 

6. 22 ±0. 12 

4.06 

9.56 

8.89 

8.45 

7. 73±0. 14 

8. 55±0. 14 

8. 64±0. 14 

4.06 

9 56 

8.86 

8.41 

7. 97±0. 14 

8. 80±0. 14 

8. 89±0. 14 

5.11 

12.90 

11.84 

11.10 

10. 19±0. 22 

11. 49±0. 22 

11. 69±0. 22 

6. 10 

16.42 

15.05 

13.93 

13. 74±0. 37 

15. 68±0. 36 

16. 02±0. 37 

6.18 

16.71 

15.22 

14. 08 

13. 66±0. 35 

15. 40±0. 35 

15. 61±0. 35 

6.61 

18.37 

16.78 

15.41 

13. 91±0. 32 

15. 98±0. 31 

16. 33±0. 32 

6.65 

18.60 

17.09 

15.70 

13. 94±0. 35 

16. 17±0.34 

16. 70±0. 35 

7.62 

22.72 

20.80 

18.80 

16. 96±0. 47 

19. 70±0. 46 

20. 27±0. 47 

8.18 

25.44 

22.91 

20.52 

19. 39±0. 60 

22. 75±0. 58 

23. 34±0. 60 

8.22 

25.66 

23.20 

20. 76 

19. 25±0. 57 

22. 29±0. 55 

22. 79±0. 57 

8.33 

26.20 

23.86 

21.29 

19. 96±0. 68 

23. 59±0. 65 

24. 50±0. 68 

9.16 

30.54 

27.49 

24.16 

22. 34±0. 82 

25. 48±0. 78 

26. 28±0. 82 

9.20 

30.77 

27.84 

24. 43 

23. 05±0. 87 

26. 35±0. 83 

27. 20±0. 87 

9.37 

31.76 

28.85 

25.23 

23. 14±0. 64 

27. 18±0. 62 

27. 91±0. 64 

10.18 

36.96 

33.14 

28. 46 

25. 83±0. 92 

31. 01±0. 89 

32. 12 ±0. 92 

10.81 

41.45 

37.56 

31.73 

32. 11±1. 5 

37. 15±1. 4 

38. 75±1. 5 

10. 98 

42.73 

38.02 

32.00 

27. 83 ±1. 5 

32. 30±1. 4 

33. 79±1. 5 

10.98 

42.73 

37.94 

31.93 

31. 07 ±0. 99 

35. 26±0. 95 

36. 11±0. 99 

11.77 

49.36 

43.34 

35.66 

34. 22±1. 5 

39. 38±1. 4 

40. 62±1. 5 

12.22 

53.44 

47.76 

38.69 

37. 25±1. 3 

43. 51±1. 2 

45. 37±1. 3 

12.58 

56.87 

50.25 

40.22 

41. 54±1. 4 

48. 34±1. 3 

50. 00±1. 4 

12.74 

58.41 

52.85 

42.17 

41. 50±2. 2 

49. 54±2. 1 

52. 64±2. 2 

12.88 

60.07 

52.89 

41.91 

43. 09±2. 4 

50. 11±2. 3 

51. 84±2. 4 


Experimental delayed critical height, 26. 94 cm; calculated delayed critical height, 26. 80 cm. 
^Experimental errors based on counting statistics only. 
c He 3 detector data, all other data with BF g detector. 
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quite apparent thatlPmuch greater scatter of data occurs for all three methods than re- 
sults from the inhour method. This scatter appears to be more a function of the delayed- 
neutron "area” than the prompt-neutron "area. " In particular, the cases with greater 
shutdown worths suffer from poor statistics. Determination of the channels over which 
the background is constant is also an important factor. This determination is done by the 
computer code GRIPE II, as discussed in Analysis of Experimental Data, but judgment 
on the part of the experimenter is necessary to assure proper use of the code. 


CONCLUSIONS 

A number of conclusions result from this study of an unreflected homogeneous, ther- 
mal reactor from delayed critical to about $50 shutdown: 

1. Multigroup-transport- theory- calculated values of the fundamental prompt- mode 
decay constant agree with experimentally determined values to within 5 percent from de- 
layed critical to $5 shutdown; from $5 to $50 the calculated values exceed experiments 
progressively from 5 to 15 percent. 

2. All methods for determining reactivity experimentally agree to about ±5 percent 
of the average of all methods for configurations up to $6 shutdown. 

3. Experimentally determined reactivities based on the inhour method, which uses 
the experimental value for the fundamental prompt-mode decay constant and a calculated 
parameter, had the best agreement with calculations over the range of measurement. 

4. The method of Simmons and King gave experimental values for the reactivity 
which differed considerably from the calculated values over most of the range of solution 
heights considered. This method fails because the value of the ratio of the effective 
delayed- neutron fraction to the generation time is not constant with increasing subcriti- 
cality. 

5. The Garelis-Russell and Gozani area-ratio methods predicted experimental values 
for the reactivity in agreement with each other and in reasonable agreement with calcula- 
tions. 

6. The area-ratio method of Sjostrand for determining experimental reactivity was 
not in agreement with calculated values over the range of the measurements . This 
method failed because of its more approximate basis as compared to the other two area- 
ratio methods. 

7. All of the area-ratio methods for determining reactivity are sensitive to an accu- 
rate experimental determination of the equilibrium delayed-neutron background. 

Thus, of the five methods considered for obtaining the shutdown reactivity of the sys- 
tem studied, the inhour method and the Garelis-Russell and Gozani area-ratio methods 
gave comparable results and agreed reasonably well with calculations. However, all the 
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area-ratio methods were sensitive to an accurate determination of^Rfe delayed-neutron 
background and exhibited some scatter in the data over the range of measurements con- 
sidered. In addition, the area-ratio methods are subject to kinetic distortion factors 
which are quite difficult to compute in practice. The determination of shutdown reac- 
tivity by using the inhour method appears to be the most satisfactory over the range of 
measurements considered. This method does require a calculated parameter, namely 
the ratio of the effective delayed-neutron fraction to the generation time, which, however, 
is relatively easy to obtain. The use of the inhour method has another significant advan- 
tage over the various area-ratio methods in that the delayed-neutron background is not 
required, thus considerably simplifying and shortening the run time for an experiment. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, August 14, 1969, 

129-02. 
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APPENDIX A 
SYMBOLS 


V> 

b n 

c 4 (r,t) 


dE 

dr 


dt 

df? 

E 

F k 



F . • 
nj,i 

F oo 

fj(E)dE 


f k (E) 

f p (E)dE 

f s (E) 

G(E* - E)dE 


H(E' - E)dE 

K eff 

Z k 

Z 0 


j n time -dependent coefficient in expansion of (F, E,£2, t) (eq. (B31a)) 

hr 

coefficient in expansion of time integral of delayed- neutron flux (eq. (B39)) 
i L delayed- neutron precursor concentration per unit volume at position 

t n 

r at time t, no. /(cm°)(sec) 
differential energy element, eV 

3 

differential volume element, cm 
differential time element, sec 
differential solid angle, sr 
neutron energy, eV 

integral over reactor volume of product of Pg(F) and P^(F) (eq. (B16a)) 


integral over reactor volume of product of P* j(F) and P^(F) (eq. (B16b)) 

integral over reactor volume of product of P^ + (F) and Pj?(F) (eq. (B45a)) 

integral over reactor volume of product of P® + . ff") and Pf(F) (eq. (B45b)) 

n, x j 

integral over reactor volume of Pq + (F) and Pq(F) 

fraction of neutrons belonging to delayed group i in energy range dE 
about energy E 

kinetic fission spectrum (eq. (B7)) 

fraction of prompt neutrons born in a fission process occurring at any 
energy and which are found in the energy range dE about energy E 

static fission spectrum (eq. (B8)) 

probability that a neutron at energy E' will, after undergoing an (n, 2n) 
reaction, be in the energy range dE about energy E 

probability that a neutron at energy E* will, after an inelastic scattering 
collision, be in the energy range dE about energy E 

neutron multiplication factor 

neutron lifetime for k in prompt mode, sec 

neutron lifetime for fundamental prompt mode, sec 
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N(a k ) 

N d (r) 

N d r,t) 

N?(r) 

K < ? > 

N p ^) 

N p r,t) 

N§(r) 

N§(r,t) 


number of delayed-neutron emitters for a given fissionable isotope 

normalization factor for kinetic fission spectrum (eq. (B9)) 

total number of delayed neutrons per unit volume at position r, 
no./cm^ (eq. (B70b)) 

total number of delayed neutrons per unit volume at position r" 

O 

and at time t, no. /(cm ) (sec) (eq. (B50a)) 

"Hi 

total number of neutrons of j 1 prompt mode per unit volume at 
position r, no. /cm (eq. (B56b)) 

total number of neutrons of n th static mode per unit volume at po- 
sition r , no. /cm^ 

total number of prompt neutrons per unit volume at position r, 
no./cm^ (eq. (B70a)) 

total number of prompt neutrons per unit volume at position r, 

Q 

and at time t, no. /(cm ) (sec) (eq. (B56a)) 

total niunber of neutrons of fundamental prompt mode per unit vol- 
ume at position r, no. /cm 0 (eq. (B52a)) 

total number of fundamental prompt-mode neutrons per unit volume 
at position r and at time t, no. /(cm ) (sec) (eq. (B74)) 


n(r, E, fT, t)dF dE d£T probable number of neutrons belonging to space volume element 

dF about position r , having energies in dE about energy E, 
and being in the element of solid angle dS? about at time t 

P(E’ — E;jUj j )dE probability that a neutron at energy E' will, after an elastic 

scattering collision, be in the energy range dE about energy E 
and with the scattering angle being through angle cosine in 
the laboratory coordinate system 

P d (F) number of fission neutrons of all energies produced per unit vol- 

ume at position r for delayed-neutron part of directional neu- 

o 

tron flux density, no. /(cm ) (sec) (eq. (B64a)) 

P?(F) number of fission neutrons of all energies produced per unit vol- 

ume at position r for j** 1 prompt mode (eq. (B44a)) 

P. (r) number of fission neutrons of all energies produced per unit vol- 

K 4-1-. o 

ume at position r fork prompt mode, no./(cm )(sec) 

(eq. (B16c)) 

P S (F) number of fission neutrons of all energies produced per unit vol- 

n ^ 4-i_ 

ume at position r for n static mode (eq. (B44b)) 
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^1111 

pf(F) 

fission-spectrum-weighted importance function at position F for 
the n** 1 static mode (eq. (B44c)) 



i tn delayed-neutron-spectrum-weighted importance function at po- 

j-i* 

sition F for n static mode (eq. (B44d)) 


p p ^) 

number of fission neutrons of all energies produced per unit vol- 
ume at position F for prompt- neutron part of directional neutron 
flux density, no. /(cm^)(sec) (eq. (B64b)) 



static-fission-spectrum-weighted importance at position F 
(eq. (B16d)) 


K,i & 

i LI delayed-neutron-spectrum-weighted importance function at po- 
sition F (eq. (B16e)) 


Q(F, E, t)dF dE dfS 

number of neutrons produced by a fixed source in space volume 
element dF about position F, having energies in dE about 
energy E, and being in the element of solid angle df? about 
at time t 


Q. 

j 

coefficient defined by eq. (B38) 


(^(F, E, f?)dF dE dj? 

AT_ 

neutrons produced by a fixed source for k time mode belonging 
to space volume element dF about position F, and with energies 
in dE about energy E, and in the element of solid angle dfl 
about 0 per sec (eq. (B5b)) 


Q 0 

coefficient for fundamental mode defined by equation (B38) 


R 

repetition rate for pulsed-neutron source, sec - * 


n 

position on surface of reactor system 


r 

position vector 


S 3 

j** 1 coefficient in expansion of Q(F, E, JT,t) (eq. (B31b)) 


T d (t) 

time -dependent function used in separating time-dependence of 
delayed-neutron directional flux density from its space-energy- 
angle dependence (eq. (B60b)) 


T p (t) 

time -dependent function used in separating time dependence of 
prompt- neutron directional flux density from ts space-energy- 
angle dependence (eq. (B60a)) 


t 

time, sec 


v(E) 

neutron speed corresponding to energy E, cm/sec 
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X 


z 00 ^) 

z£°(r) 

Z GR (r) 


z£ R (r) 

Z S0 (r) 



(?) 


a? 


a 


a- 


a 



P 


0 

h 

% 


parameter of Gar elis -Russell area- ratio method, defined as a solution of 
eq. (B54) 

spatially dependent kinetic distortion factor for extrapolated area-ratio method 
of Gozani (eq. (B52b)) 

spatially dependent kinetic distortion factor for Gozani method for fundamental 
prompt and delayed modes only (eq. (B53)) 

spatially dependent kinetic distortion factor for Garelis-Russell method 
(eq. (B58)) 

spatially dependent kinetic distortion factor for Garelis-Russell method for 
fundamental prompt and delayed modes only (eq. (B59)) 

spatially dependent kinetic distortion factor for area- ratio method of Sjostrand 
(eq. (B72)) 

spatially dependent kinetic distortion factor for Sjostrand method for funda- 
mental prompt and delayed modes only (eq. (B72*)) 

f Vi — 1 

decay constant for j prompt mode, sec" 

decay constant for k Ln mode, sec" 

tli -i 

decay constant for k delayed mode, sec" 

decay constant of fundamental prompt mode, sec"'*' 

experimentally measured value of fundamental prompt-mode decay constant, 
sec"* 

experimentally measured value of fundamental prompt-mode decay constant 
for reactor system at delayed critical, sec"* 

total delayed-neutron fraction per fission for given fissionable 
m 

isotope, p Pj 
i=l 

effective delayed-neutron fraction for total prompt flux (eq. (B68)) 

J--U 

fraction of neutrons emitted per fission and belonging to i delayed group for 
given fissionable isotope 

■fli -j-u 

effective delayed-neutron fraction for i Ln delayed group and k prompt mode 

XT. 

total effective delayed-neutron fraction for k u prompt mode (eq. (B19)) 

xu xu 

total effective delayed-neutron fraction for n n delayed mode and j n prompt 
mode (eq. (B47)) 


28 



*nj,i 

* 

^0 

^00 

A k 

A 0 

X i 

m l 

v(E) 

P s 

P s ($) 

P^($) 

Ps R ($) 

P^°($) 

P^($) 

P® K ($) 

pg°($) 

Ps 

S a (E) 

S es( E ) 

S f (E) 

S in (E) 


ii* 

effective delayed- neutron fraction for i delayed-neutron group corresponding 
to n** 1 delayed mode and j** 1 prompt mode (eq. (B46)) 

fraction of prompt neutrons emitted per fission for given fissionable isotope 

total effective delayed-neutron fraction for fundamental prompt mode 

total effective delayed-neutron fraction for fundamental delayed and prompt 
modes, 0 OO = 0 O 

generation time for k 1 prompt mode, sec 
generation time for fundamental prompt mode, sec 

"Hi — "I 

decay constant for i ul group delayed-neutron precursor, sec" 

cosine of angle in laboratory coordinate system between initial and final direc- 
tion of motion of elastically scattered neutron 

number of neutrons produced per fission at energy E 

static reactivity 

static reactivity in dollars 

calculated reactivity in dollars based on either equation (B20a) or (B20b) 

experimentally determined reactivity in dollars based on Garelis-Russell 
method (eq. (B77)) 

experimentally determined reactivity in dollars based on Gozani method 
(eq. (B76)) 

experimentally determined reactivity in dollars based on inhour method 
(eq. (B75)) 

reactivity defined for n in and m Ln static modes 

experimentally determined reactivity in dollars based on Simmons-King 
method (eq. (B79)) 

experimentally determined reactivity in dollars based onSjostrand method (eq. (B78)) 
reactivity for fundamental delayed mode 

macroscopic neutron absorption cross section at energy E, cm"'* 
macroscopic neutron elastic scattering cross section at energy E, cm * 
macroscopic neutron fission cross section at energy E, cm"* 
macroscopic neutron inelastic scattering cross section at energy E, cm 
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s n2n< E > 

S t (E) 

S S (E' - E; m l ) 

T 

0(r, E,f2,t) 

0 d (r,E,f2, t) 
0 p (r,E,f2,t) 
</> d (r,E,f2) 

<^(r, E) 
<P k (r,E) 
<P k (r,E,fS) 

<p£(r,E,sl) 

</>g(r,E, f!) 

.pP+^E,^) 

<P„(r,E) 
<p®(r,E, jf) 


macroscopic neutron (n, 2n) reaction cross section at energy E, cm” 

macroscopic neutron total cross section at energy E, cm - *, 

S t (E) = S a (E) + S eg (E) + S in (E) + S n 2 n (E) 

generalized scattering kernel (eq. (BIO)) 

width of a source pulse, sec 

directional flux density at point r and time t for neutrons of energy 

— *• O 

E traveling in direction SI, no/(cm ) (sec) (unit energy) (sr), 

0(r, E, f?, t) s v(E)n(F, E, S?, t) 

delayed component of 4>(r, E,f2, t), no. / (cm 2 ) (sec) (unit energy )(sr) 

prompt component of <p(r, E, fi, t), no. / (cm ) (sec) (unit energy)(sr) 

directional delayed-neutron flux density at point r for neutrons of 
energy E traveling in direction O, no./(cm ) (sec) (unit energy) (sr) 
(eq. (B60b)) 

scalar neutron flux for j n prompt mode at position r for neutrons of 

p 

energy E, no. /(cm ) (sec) (unit energy) 
scalar neutron flux for k u 1 time mode at position r for neutrons of 

o 

energy E, no. /(cm ) (sec) (unit energy) 
directional flux density at point r for neutrons of energy E traveling 

ru O 

indirection S2 for k time mode, no. /(cm ) (sec) (unit energy) (sr) 
(eq. (B5a)) 

directional flux density at point r for neutrons of energy E traveling 

n i-U 

in direction it for k n delayed time mode, 

p 

no. /(cm ) (sec) (unit energy)(sr) 

directional flux density at point r” for neutrons of energy E traveling 

t fVi 0 

in direction for k prompt time mode, 

9 

no. /(cm ) (sec) (unit energy) (sr) 

j.’l 

adjoint flux corresponding to k n prompt time mode and which gives the 
importance of neutrons at position r which have energy E and are 
going in direction 

scalar neutron flux for n** 1 static mode at position r for neutrons of 

o 

energy E, no. /(cm ) (sec) (unit energy) 

j_i_ 

directional n n mode static flux density at point r for neutrons of 
energy E traveling in direction and corresponding to the reac- 
tivity no. /(cm 2 ) (sec) (unit energy) (sr) (eq. (B12)) 
s 
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(pf. + (r, E) scalar adjoint neutron flux giving importance of neutrons at position r and 
having energy E for the n static mode 

® + (r, E, to) static n th mode adjoint flux defined as a solution to equation (B13) and which 
gives the importance of neutrons at position r which have energy E and 
are going in direction to 


<p p (r,E, f!) 


directional prompt-neutron flux density at point r for neutrons of energy 
E traveling in direction no. /(cm 2 ) (sec) (unit energy) (sr) (eq. (B60a)) 


<P s (r,E,S2) 


<(r,E) 

<?g(r,E, f?) 


to 


directional flux density at point r for neutrons of energy E traveling in 
direction f! corresponding to the reactivity p_, 

o b 

no. /(cm ) (sec) (unit energy) (sr) (eq. (B12)) 

static adjoint scalar flux giving the importance of neutrons at position r 
and having energy E 

static adjoint flux defined as the solution to equation (B13) and which gives 
the importance of neutrons at position r which have energy E and are 
going in direction S5 

unit vector in direction of motion of neutron travel 


to + 

to~ 

$ 


directions which are outwardly directed on surface of reactor 

directions which are inwardly directed on surface of reactor 

unit of reactivity defined such that $1 represents the amount of reactivity 
required to change a delayed critical reactor to prompt critical 
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APPENDIX B 


THEORETICAL ANALYSIS 
Conventional Static Reactivity 

The theory underlying pulsed-neutron-source experiments is based primarily on the 
a na lysis given in references 1 to 9 and 12. This appendix is based on the review article 
of Preskitt et al. (ref. 4) and gives the various equations in conventional rather than op- 
erator form. Also the working forms of the various equations are indicated. The equa- 
tions for the area-ratio method of Sjostrand are based on reference 12. 

To provide a theoretical basis for understanding pulsed- source neutron experiments, 
consider the time-dependent transport equation with time -independent cross sections and 
a single fissionable isotope. With these restrictions, the number density of neutrons at 
position r having direction H and energy E at time t n(r,E,fi,t) is given by 


9n(r,E,f?,t) + ^ . Vv(E)n(r, E, fj, t) + S t (r, E)v(E)n(r, E, f?, t) 

at 

= Q(r,E,^,t)+ JJ S eg (r,E , )v(E')n(r,E , ,j5 , ,t)P(E' - E;p L )dE' dJ5’ 

+ 1- f f S.^E'ME’jnfcELfSLOHCE’ - E)dE* dST 
An JJ in 

+ J JJ S n2n (r,E , )v(E')n(r,E , ,fJ',t)G(E’ - E)dE* dn’ 


f„( E ) rr 

/3 p ^-— JJ i'(r,E')S f (r,E>(E')n(r,E\ j?',t)dE’ dsT 
m , 

Z f. (E) A (t-t T ) rr 

- — I dt'e / / !^(r, E f )Sj(r, E , )v(E , )n(r, E',^' ,t')dE' df5' 

i=l Cun 


The notation used for equation (Bl) is standard, and all symbols are listed in appendix A. 
The last term on the right side of equation (Bl) is given by Gross and Marable (ref. 19). 
In equation (Bl) the following normalizations are required: 
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/ f (E)dE = 1.0 
* / E p 


(B2a) 



(i = 1, . . . , m) 


(B2b) 


f dti f P(E f — E; jn T )dE = 1.0 
•/j! /E 

r e ’ 

H(E' — E)dE = 1.0 

r E’ 

^ G(E' — E)dE =2.0 


(B2c) 

(B2d) 

(B2e) 


The total delayed- neutron fraction (3 for a given fissionable nuclide is given by 

m 

(B3a) 

i=l 


and 

/3 +0 = 1.0 (B3b) 

Define the directional flux 0 in equation (Bl) in terms of the product vn, that is, 
let 

0(r, E, fj, t) = v(E)n(r, E, t) (B4) 

Now solutions to equation (Bl) are sought which vary exponentially in time, namely, 

0(r, E, £?, t) = <P k (r, E, S?)exp(a! k t) (B5a) 

Q(r, E, t) = Q k (r, E, f?)exp(o! k t) (B5b) 

The expressions given by equations (B4), (B5a), and (B5b) are substituted into equa- 
tion (Bl) and give 


33 



f2 • ^<P k (r, E, 55 ) + js t (r, E) + <? k (r, E, S?) 

= Q^.(?;E,Sf)+ JJ SgfcE' -E^L^k^ES^JdE'dS' 


fk(E) 
47 r 


/y* ^(r, E*)S f (r, EO^jjff’, E'jfJ^dE’ dfT (B6) 


Here f k (E) defines a kinetic fission spectrum given by 

m m 

v e > = w e) + 


y hViW f (E) y a k¥i< E > 

Zj “k +x i s Zy “k +x i 

i=l i=l 


(B7) 


where f_(E) is the static fission spectrum given by 

S 

m 


f s (E) = /3 1(E) + ^ ^fj(E) 
i=l 


(B8) 


In the limit as « k approaches zero, the spectrum f k (E) approaches the steady-state 
fission spectrum f„(E). Since both the prompt and delayed fission spectra are normal- 

b 

ized to 1, the kinetic fission spectrum f k (E) is normalized to 


m 

f f k (E)dE = 1 - S' • — - - s N(Q! k ) 

J E k Z^“ k + X i 


(B9) 


i=l 


In equation (B6) the generalized scattering kernel has been defined as 


S s (r, E* - E; m l ) = £ e s (r > E ’ )P(E ' E > ^l) + s in (r ’ E,) E) + S n2n^ E>) ^ ~ E) 

rxTT Tt 7i 

(BIO) 


Equation (B6) is the expression satisfied by all exponentially varying solutions to 
equation (Bl) and contains no references to the reactivity of the system. Reactivity 
must then be defined in terms of modal quantities, such as « k and <? k , which can be 
measured. The reactivity used here is the static reactivity p g , which is obtained in 
conventional reactor calculations. That is, let 
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I 


= 


K eff - 1 


K 


eff 


(Bll) 


where is the multiplication factor . Thus, the static reactivity p g is obtained as 

a solution to the equation 


S2- ^<P s (r, E,f?) + S t (r, E)<p g (F, E,S?) 

= JJ S g (F, E' — E; n-jjcp s (p, E f , n'jdE* dfT 


+ (1 - pJ 


f s ( E > 

477 


JJ v(?,E’)X f (?,E')<p s (?,E\$')dE' dif’ (B12) 


Equations (B6) and (B12) are similar in form, with equation (B12) being the equation 
solved for various reactor analysis problems. However, solutions to equation (B6) may 
be obtained by using equation (B12) if a term a^/v{E) is added to E^.(F, E) and if the 
spectrum (1 - p g )f g (E) is replaced by the spectrum ^(E). The source term Q^(F, E,S?) 
in equation (B6) can be added to equation (B12) if needed. 


Inhour Equation 


To relate the static reactivity p g to the modal quantities and cp-^ write the ad- 
joint of equation (B12), that is, 

-ff • ^<p g (r, E, if) + S t (r, E)^ g (r, E,S?) 

* JJ S g (r, E - E’;M L )^ g (r, E',f?')dE' dS?’ 

v(jc, E)2-(F, E) rr 

f - JJ yEV^E’^’ME’ dti' (B13) 


+ ( 1 -P S ) 


47T 


Now multiply equation (B6) by <P+{r, E, 0) and integrate over all dF dE dfS; then multiply 
equation (B13) by ^^(F, E,S?) and integrate over all dF dE d?S; now subtract the second 
equation from the first. The result will be the inhour equation: 



(B14) 
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where the various quantities of equation (B14) are defined as 

A k = — fffaK E, SS) — — <Py.(r, E, n)dr" dE df5 
Ft, J J J v(E) 


(B15a) 




(B15b) 


along with 


F k = / 

(B16a) 

F k,i = / 

(B16b) 

(r) = J V<p, E)S f (r, E)^ k (r,E)dE 

(B16c) 

Pg(r) = y*f s (E)<?7 + (F,E)dE 

(B16d) 

Pg >i (r)= / f i (E)^ + (r,E)dE 

(B16e) 

<P k (r,E) <? k (r,E, j?)df? 

(B16f) 

«pg(r,E) = f <Pg(r,E,f£)dn 

(B16g) 


To obtain equation (B14), the fixed source term in equation (B6) has been set to zero. 
The terms involving ^<? k (r, E, £?) and ^<pg(r, E, ST) have been eliminated by using 
Green's theorem along with the boundary conditions 


(p t_ (R, E , j? ) — 0 


<^(R,E,?T) = o 


(B17a) 

(B17b) 


with R being on the surface of the reactor system. The direction refers to all in- 
wardly directed directions on the reactor surface, while refers to all outwardly di- 
rected directions. 

If is the decay constant of the fundamental prompt mode, equation (B14) becomes 


p s - q-Pa 0 + 


(I «g| > i = l,m) 


(B18) 
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The quantity /3 Q is the effective delayed fraction for the fundamental prompt mode. For 
the k th prompt mode, the effective delayed- neutron fraction is given by 

m 
i=l 

Similarly, if A k is defined as the generation time for the k prompt mode, 
generation time for the fundamental prompt mode. 

In terms of dollars, the calculated reactivity (eq. (Bll)) is given by 

(B20a) 

*0 K ef#0 

The reactivity in dollars as calculated by using equation (B18) is 


(B19) 

Aq is the 


P°($) =^=^-+ 1 

^0 ^0 

A o 


th 


If the neutron lifetime, Z k for the k prompt mode is defined as 


l k ~ K eff A k 


(B20b) 


(B21) 


equations (B14), (B18), and (B20b) may be rewritten in terms of the neutron lifetime, 
namely, 

/ 

m 


p s = 


«kf— + y 

\ K eff Z-/*k + X i 

\ i=l 


(B14 1 ) 


„ “<V T 

p a = — + 00 


K 


eff 


(B18’) 


P Q ($) =' 


a*. 




+ 1 


(B20b') 


eff p 0 


Equation (B18) shows that calculated values of A fl and /3 Q are required to obtain 
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the reactivity p g from experimental values of the fundamental prompt- mode decay con- 
stant Q!q. Two static calculations are required so that equation (B18) may be compared 
to experiment. The first static calculation yields a solution for equation (B12) and gives 
p g . The second static calculation gives a solution to equation (B6) for with the fixed 
source term set to zero. The spectrum f k (E) for this calculation, given by equation (B7), 
goes to 


f k ( £ ) - |8 p fpCE) (|ajj}| » i = 1, m) 


(B22) 


along with 


1 - P s = /3 p (B23) 

Then, when calculated values of afi and p g are known and Bq is obtained by some other 
calculation, equation (B18) can be solved for Aq. Finally, when experimental values of 
Q!q and the calculated values of Aq and Bq are used, "experimental* ' values of the 
reactivity p g can be obtained by using equation (B18). 

Another procedure for obtaining all of the parameters of equation (B18) makes use of 
the solution to equation (B6) as previously described and the solution to the static adjoint 
equation given by equation (B 13). Then Aq is determined from equation (B15a) and /3 q 
from equations (B15b) and (B19). 


Reactivity from Modal Areas 


Expressions can be developed to relate the static reactivity to the modal fluxes and 
the decay constants. This approach is called a modified pulsed-source technique and the 
best known methods are those of Gozani (ref. 2), Garelis and Russell (ref. 3), and 
Sjostrand (ref. 1). 

The distinction between prompt and delayed modes can be seen from equations (B6) 
and (B7). For prompt modes, |a! k | » for all delayed-neutron groups, and the source 
spectrum is given by equation (B22). The equation for the prompt modes is then 


S2 • V<pP(T, E, S2) + 


S t (r, E) + 


af 


v(E) 


pg(r,E,s5) 


= JJ E s (r, E' - E; p L )^g(r, E*,f2 , )dE* dfT 

b f (e) r r 

+ -2-2 // v(p, E')S f (r, E')<pP(r, E',fr)dE' dS2' 

4ttK .. JJ 1 K 


(B24) 
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For delayed modes, the decay constants are known to cluster about the precursor 
decay constants Xp and the term o;^/v(E) in equation (B6) is negligible in comparison 
to the remaining terms on the left side of the equation. Thus, the delayed modes are 
given by 


& ’ > E, ^2) + 2j.(r, E, 

= ff 'L b &,V' ~ EjML^fcE’.sS'JdE* diT 


fv(E) 


4ttK 


eff 


j"J v(jr, E’)S f (r, E’)^^?, E', fT)dE' dH' 


(B25) 


Comparing equation (B25) with the static equation (eq. (B12)), it can be seen that the 
solutions for the delayed modes differ from the static flux only in the effect of a different 
source neutron energy spectrum. The prompt modes differ from these delayed modes by 
the additional influence of the term cn j^/v (E ) on the neutron spatial and energy distribu- 
tion. This difference between the prompt and delayed modes is referred to as kinetic 
distortion. 

The small difference between the static and delayed modes caused by the different 
energy spectrum of source neutrons is ignored. The prompt modes are thus described 
by equation (B24) and the delayed modes by equation (B12). By writing the adjoint equa- 
tions corresponding to equations (B12) and (B24), it can be shown that the following or- 
thogonality relations hold: 


He - ®?) fff ^ dE <^ + (r>E, ??)-!_ <pP(r,E, 


m 


= 0 


(B26) 


n m 

p s -p s 

An 


/d? JJ dE dfS <p^ + (r,E,ft)f s (E) 


JJ i^E'JSjfcE'^fcE’^'ME' d^’ =0 


(B27) 


Equation (B26) states that the prompt adjoint flux is orthogonal to the prompt-neutron 
density for different prompt modes. Equation (B27) states that the delayed adjoint flux 
is orthogonal to the delayed fission source for different delayed modes. 

A pulsed- neutron transient is now described by assuming that the prompt modes are 
excited directly by the external source pulse and that the delayed modes are excited by 
the decay of precursors formed throughout the transient. This is the approach followed 
in reference 7 and extended by Preskitt et al. in reference 4 to display the role of higher 
exponential modes, both prompt and delayed, as well as of the kinetic distortion. For 
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this analysis, it will be convenient to rewrite equation (Bl) as (m + 1) time -dependent 
equations: 


_j_ d ± (r, E, g, t) + ^ . v0(r,E,n,t) + S t (r,E)<^(r,E,n,t) 
v(E) 3t 


= Q(r,E,&,t) + JJ S s (r,E' - E; M L )0(r, E , ,S?',t)dE’ dST 
a i (e) r r 

+ — E-B JJ v (p } E T )2 f (r, E’)^?, E', fT,t)dE' dfT 


X i C i (F,t)f i (E) 

+ (B28a) 

4n 


3C- (r, t) r r 

= i8. y J v(r, E')S f (F, E')$(F, E', f2 f ,t)dE' dS?' - A-C^t) 


(i = 1,2, . . . ,m) 


(B28b) 


Equations (B28a) and (B28b) are now written in terms of the prompt component of the an- 
gular flux <f>p(r”, E, f?, t) and the delayed component of the angular flux 0 d (T, E , t): 

i E, fi, t) _ 

— 2 + * V0 (F,E,0,t) + 2 t (r,E)<£ (r,E,S2,t) 

v (e) at p 1 p 

= Q(F,E,f!,t) + // S s (F,E' -EjqL^pfcE’^tME'dfT 
/3 f (E) /*/* 

+ -£-2 // i^E^S^E’)^ (^E'^'^dE' dfT (B29a) 

4t\ J J 1 p 


j 00 d (r, E,St,t) 


+ n • V0 d (r, E,S2,t) + S t (r,E)0 d (r,E,S2,t) 


v(e) at 

- ff 2 s (F,E’ -E; f i L )0 d (r,E , ,^ , ,t)dE’ dfT 

W«CE) 


P P V 


If.,-* 


A-C. (r, t)f- (E) 

)2 f (r, E , )0 d (r,E , ,fi’,t)dE' dfT + — - 1 -- (B29b) 


4 ir 
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(i = 1,2, . . . , m) 


(B29c) 


0C. (r,t) cr — — r — — 

— ^ = /3 i J J ^( r , E, )S f (r,E')[^ p (r,E , ,S2',t) 

+ 0 d (r,E , ,fS',t)J - XjCjC^t) 


<t>( r,E,^,t) = 0 p ( r,E,jS,t) + <j> d ( r,E,S?,t) 


(B29d) 


Now let a single source pulse occur at time t = 0 and integrate equations (B29b) and 
(B29c) from t = 0 to t = 00 . Zero initial and final conditions are assumed. If the re- 
sulting two equations are now combined by eliminating the term involving the precursor 
concentration, there results 


["■'•*«] r 0 d (r, E, J?,t)dt 

| JJ dE' dS' E g (r,E’ -* E;p L ) + JJ dE* dS?’ v&, E’)S f (r, E f ) 


•f 0 d (F,E\j5',t)dt + 


m 


^(E) 


47T 


jy dE T dst’ v(r, E')S f (r,E') 


L i=l 


f <M r,E\fT,t)dt = 0 
y 0 p 


(B30) 


This result is obtained for the case of a single pulse; however, the analysis also holds 
for a repetitively pulsed case (ref. 20). Equation (B30) represents the starting point for 
obtaining an expression where the modal fluxes are related to the reactivity for the mod- 
ified pulsed- source technique. 

Assume a square source pulse of width r, expand the prompt flux in terms of the 
eigenfunctions ^(r’, E, Vi), and also expand the external source in terms of the prompt- 
neutron densities v(E)' " 1< Pj 3 (r, E, j?): 


0 p (r, E,f?,t) = £ A.(t)<pP(r, E,fJ) 
j 


Q(r, E, Vi, t<r) = s, -i— <p?{ r, E, f?) 

3 v(E) 3 

j 


(B31a) 

(B31b) 


Substitute equations (B31a) and (B31b) into equation (B29a) and simplify the resulting 
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II 


1 


equation by using equation (B24). Now multiply the equation by cp^ + (r, E, S?) and inte- 
grate over all dF dE dS5. Finally, making use of the orthogonality relation given by 
equation (B26) yields 


///*- 


dE dO 


•I v(E) 3 


|dA.(t) „ , 

— «fA.(t) - s.f = 0 

dt 33 3 I 


(B32) 


Also, equation (B26) implies that the term in the braces of equation (B32) must be zero, 
namely, 


dAj (t) 

dt 


“jV>-s r O 


The solution to the homogeneous form of equation (B33) is 


(B33) 


A j(t) = Qj exp(Q'Pt) (t > r) 


(B34) 


where Qj is, as yet, an undetermined coefficient. 

To obtain a solution of equation (B33) for the small time interval from t = 0 to 
t = r, the following artifice is used. Convert the differential equation to an algebraic 
equation, that is, 


t 


“SV* - s i = 0 


ft = r) 


(B35) 


which gives 


A j( t) 



t-aPt 


(t ^ r) 


(B36) 


If the pulse width r is assumed to be small compared to l/\aP\, or equivalently 
|aj>t|«|, equation (B36) can be written as 


A.(t) = Sjt exp(o!?t) (t ^ r) 


(B37) 


When A j(t) from equation (B34) is used, equation (B37) evaluated at t = r gives 

Qj = TSj (B38) 

For a square source pulse, s, is the magnitude independent of t, and Q- represents the 

+K •» J 

source strength for the j prompt time mode. 

In making the expansions given by equations (B31a) and (B31b), it is assumed that the 
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prompt modes form a complete set. Thus, the small time interval during which the 
pulsed- source neutrons are slowed down to fission neutron energies is ignored (ref. 9). 
Next expand the delayed flux integral in terms of the static eigenfunctions, that is, 

f 0 d (f,E,fi,t)dt = ^b n <?p®(r,E,S2) (B39) 

® n 


Substitute equation (B39) into equation (B30). Multiply the resulting equation by 
<p^”(r, E, S2) and integrate over all dr dE dO. Making use of equation (B12) and the or- 
thogonality relation given by equation (B27) results in 


fff dF dE dn <?® + (F, E, f?)b n • 

■/// 


4rr 


ff E .^. )dE . dfT 


dr dE df! <p® + (F, E,ff) 


m 

z^ 


(E) 
4 ir 


i=l 


ff dE' dfS’ y(F, E')S|(F, E') f 


0 p (F,E\sT,t)dt 


(B40) 


If equation (B31a) is substituted into the prompt- neutron time integral of equa- 
tion (B40), 

/ °° V — \ TS- 

<£ (F,E,j?,t)dt = -\ — 1 <p?(F,E,fS) (B41) 

u P a P 3 

3 3 


when equations (B34) and (B38) are used. Substitute equation (B41) into equation (B40) 
and make use of integrals of the type given by equations (B16f) and (B16g), then, solve 
for b n : 


)Sf(F, E')^?(F, E'JdE’ 


b n = 


Jj dr dE <p® + (F,E)^-^ j v(?, E')S f (r, E’)<p®(r, E’)dE' 


(B42) 


Now use definitions similar to those given by equations (B16c), (B16d), and (B16e) in 
equation (B42): 
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(B43) 


Z— / i=l 




/<rp^ + cf)l^(r) 


Explicitly written, these definitions are 

P?(r) = / ^E)S f ^E)^P(r,E)dE 
p ^(r") = /^(F,E)S f (r,E)^®(r,E)dE 

P n +(F) = f f s (E)<P n +(F ’ E)dE 
P®^(r) = / f i (E)< ? ® + (r,E)dE 

Define modal quantities similar to those given by equations (B16a) and (B16b): 

F nj = / P® + (r ) p P(r)dr 

F nj,i = 

Now define the effective modally dependent delayed- neutron fraction as 


nj,i 


B.F . . 
_ p i n],i 


F • 
nj 


along with 


m 


^nj X ^nj,i 
i=l 


m 


E 

i=l 


Vni.i 

F ni 


Substitute equation (B45b) in equation (B43) and then use equations (B47) and (B45a) 
to obtain the following result for b n : 


(B44a) 

(B44b) 

(B44c) 

(B44d) 

(B45a) 

(B45b) 

(B46) 


(B47) 
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(B48) 



With this solution for the delayed modes, measurable quantities may be expressed in a 
a form which yields reactivity. This is done by substituting for b in equation (B39) its 
value from equation (B48). Then this resulting expression is multiplied by v(E) - and 
integrated over all energies and all values of £2. The resulting expression is given by 




f P nV ) P n<T ^ 



(B49) 


where 


d (r,t)= / dE — - — Un <£ rf (F,E,s1,t) 

J v(E )J a 

N ®(r) = / dE J— /d£? <^®(r,E,£5) 

n J v(E )J n 


(B50a) 


(B50b) 


Rearranging equation (B49) to extract the reactivity in dollars for the fundamental 
mode yields 


P s ($) 


P 


0 

s 


^OO 


Q 0 N g(r) 


a* 


J N d (r,t)dt 
i) a 


Z G0 (r) 


where 



J d£2 <pg(r,E,£2) 


(B51) 


(B52a) 
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Equation (B51) is the fundamental relation for the interpretation of the modified pulsed- 
source method known as the Gozani or "extrapolated area-ratio” method (ref. 2). The 
quantity Z (r ) is space dependent and is called the kinetic distortion factor. 

Preskitt et al. (ref. 4) make the following observations concerning the space- 
dependent correction factor Z (r): 

(1) If kinetic distortion is present, all orders of delayed harmonics appear, even 
though only the fundamental prompt mode is excited by the pulsed source. 

(2) If no kinetic distortion is present, a given delayed mode appears only if the cor- 
responding prompt mode is excited by th*e pulsed source. 

(3) The influence of the delayed harmonics is suppressed by factors Qj/Q 0 , a^/aP, 

and p 0 /p n , and the space dependence of the correction factor tends strongly to the ratio 

of the fundamental delayed and prompt modes. 

The last consideration discussed concerning the kinetic distortion factor suggests 
G0/~* 

evaluation of Z (rj for the fundamental delayed n = 0 and prompt j = 0 modes only. 
This gives Z G0 (r ) as 


GO n5(t) / P SV)P§(r)dr 
Z«l<r) = - 

N g(r) J P® + (r)P®(r)dF 


(B53) 


which can be calculated relatively easily. 

The space dependence of other modified pulsed- source methods can also be ex- 
pressed in terms of the prompt amplitudes of equation (B34) and the delayed-neutron in- 
tegral given by equation (B49). The /3/A method of GareJis and Russell (ref. 3) is ob- 
tained by assuming the following relation between N p (F, t) and N^r, t) for some constant 
x: 


/ oo 

N p (r,t)[exp(xt) - l]dt N d (r,t)dt 


(B54) 


The equations (B31a) and (B34) are used to obtain 


N p (r,t) = ^ Qj exp(aPt)N?(r ) 
3 


(B55) 



where 


N p (r,t) = 



dfJ 0 (r, E,^,t) 

r 


(B56a) 


N?(r ) = f dE — / dfi <$(r, E, f?) (B56b) 

1 J v(E) J ^ 

Now substitute equations (B49) and (B55) into equation (B54) Perform the integra- 
tion over time, and rearrange the equation in a form such that the reactivity may be ex- 
tracted. Then, solving for the reactivity in dollars yields 



Similar remarks apply to equation (B57) concerning the spatially dependent kinetic dis- 
tortion factor Z GR (r ) as apply to equation (B51) and Z^^). Evaluating the kinetic 
distortion factor Z (r ) for the fundamental delayed n = 0 and prompt j = 0 modes 
only gives 


Z GR (r) 


N*(r) /P 5^*8^ 
N 0 (r) f P® + (r)P®(r)dr 


(B59) 


Thus, we have established that Z GR (r ) = Z^°(r ). In other words, the kinetic distortion 
correction factor based on the fundamental delayed and prompt modes only is the same 
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for the Gozani and Garelis-Russell modified pulsed-source techniques. 

The area-ratio method of Sjostrand can be derived by using equations (B29b) and 
(B29c). As before, let a single source pulse occur at time t = 0 and integrate equa- 
tions (B29b) and (B29c) from t = 0 to t = °°. Zero initial and final conditions are as- 
sumed. Combine the resulting two equations by eliminating the term involving the pre- 
cursor concentration to obtain equation (B30). Now assume that the delayed- and 
prompt-neutron fluxes are separable into a time- dependent part and a position-energy- 
angle part. That is, let 


$p (r, E , t ) = tp p (r ; E , ft )T p (t) (B60a) 

0 d (r,E,f2,t) = <p d (F,E,fS)T d (t) (B60b) 

Substitute equations (B60a) and (B60b) in (B30). This gives 


T p (t)dt n- V<P d (r ) E ) n)+ S t (F,E)< ? > d (r;E ) Sj)- f J X S <F, E'-E; ll L )V A fc E', S2’)dE' d5* - f f vp, E')<p d <?,E',fi')dE’ 

m 

Z j8.f.(E) rr _ 


f T d« dt 


Before it was assumed that the delayed modes satisfy equation (B12), and now it is 
assumed that ® d (r, E, fi) will also satisfy equation (B12). This equation is now substi- 
tuted into equation (B61) to give 


/ T (t)dt 
‘ / 0 p 

T d (t)dt 

0 



Ps ^ (E) ff v(r, E')S f (r", E’)<^ d (F, E', fT)dE' dfT 


m 

Z /3.f. (E) rr 

JJ p<r,E')E f (jr,E')<p p {jr,E',f!i')aE 1 dj? 

i-1 


(B62) 


Multiply the numerator and denominator of equation (B62) by the adjoint flux <pt(r, E,J5) 

s 

which is a solution of equation (B13), and integrate over all energies E and over all di- 
rections The result is 
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/ OO 

u v» 

f T d<‘> dt 


If f s (E)^g(r, E,fJ)dE df£ If v(p, E')2 f (r, E')^ d (F, E*, fJ')dE' dfT 


Z //“ 


(pt.(jr, E, ?X)dE dS? 

b 


II ^ 


)S f (r, E')^p(r, E’, fT)dE' dsT 


Analogous to equation (B16c) define 

P d (r ) = //,(r,E’) Sf (F,E> d ^E-,^)dE-d^ (B64a 

P p (r) = ff v(?,E')i: f (r,E')<p p (r,E\ti')dE' d& (B64b 

Now integrate equation (B63) over Then, use equations (B16d), (B16e), (B64a), and 
(B64b) to obtain 


{ T P (t>dt _ 

°o m 

/ Td(t)dt 


Integrate the right side of equation (B65) over the volume of the system: 


jf y*> dt -p a / p s (F)p d< r ) dF 
/ T d(« dt /w** - 


Then, using equations (B15b) and (B19) and assuming that P d (F) « P (r ) yields 


Pe($)=^=-fco 


-f“ V> dt / jf T p< t)dt 


— >*oo /• 

* ^ T d^^ dt / Ps^ P d^> dir J 0 T d^ dt 


I 


J 


where /3 is defined as 


0 = 


E *i ff P+ s, i(r)P p (r)dF 
i-1 

ffp*(r)P p (r)6r 


(B68) 


Multiply equations (B60a) and (B60b) by l/v(E) and then integrate over all energies 
E and all directions of motion to. The result is 


N p(r,t) = N p(D T p (t) (B69a) 

N d (r,t) =N d (r)T d (t) (B69b) 


where 


Np(r") and N^(F) are defined by the equations 


N p^) 

N d & 


/ dE -J— / dfj 

J v(E) J 

/ dE J— / dft 
J v(E) J 


<p p (r,E,£?) 
cp d (r, E,fi) 


(B70a) 

(B70b) 


Substitute equations (B69a) and (B69b) into equation (B67). The result is the area- 
ratio method of Sjostrand, that is. 


P s ($) = - 


/ oo 

y r ’ t)dt 

N (r) r° 

P J o N d (r,t)dt 


(B71) 


Since the right side of equation (B71) is a constant, the ratio N d (r )/N p (r ) is a spatially 
dependent factor which is nonunity if the delayed-neutron fluxes differ from the prompt- 
neutron fluxes. Let this ratio, the kinetic distortion factor, be defined as Z S ^(F), that 
is, 

qn N H (r) 

Z S0 (r)=— (B72 

N p (F) 

so that equation (B71) becomes for the Sjostrand method 
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(B73) 


P s ($) = - 


jT N p^ t)dt 

r N d (?-,t)dt 


Z S0 (r ) 


An approximate value for the kinetic distortion factor for the Sjostrand method may 
be obtained by using fundamental mode values of N^(F) and Np(r"). That is, let 


7 S0 

Z 0 


(?) 


Npf) 

Ng(?) 


(B72') 


Calculational and Experimental Procedure 

Since some of the experimental procedures use calculated parameters, the basic 
calculations are discussed first. The solution to the various equations can be obtained 
by using, for example, the S n method as described in reference 14. The steps in the 
calculational procedure (Cl to C4) are as follows: 

(Cl) Solve equation (B24) for the fundamental prompt-mode flux together with equa- 
tion (B13) for the fundamental static adjoint flux. These two calculations performed to- 
gether give 

(a) Fundamental prompt-mode decay constant, 

(b) Static multiplication factor, ^ 

(c) Total effective fundamental mode delayed- neutron fraction, /3 q 

(d) Neutron generation time for fundamental prompt mode, A Q 

(e) Fundamental prompt-mode flux, </?]^(r, E) 

(f) Fundamental delayed-mode adjoint flux, <pg + (r, E) 

(g) F 00 = /pg + (r)Pg(r)dr 

(C2) Solve equation (B12) for the static flux together with equation (B13) for the static 
adjoint flux. These two calculations performed together yield 

(a) Static multiplication factor, 

(b) Fundamental delayed-mode flux, <pg(r,E) 

(c) Fundamental delayed-mode adjoint flux, <pg + (r,E) 

(d) Fq Q = f Pg + (r )P®(r )dr 

(C3) From the calculations described in step (Cl), quantities can be obtained which 
may be compared to experimentally measured values. These are 
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(a) Fundamental prompt-mode decay constant, which is compared to an experi- 

mental determination, 

(b) Reactivity in dollars computed by the equation (B20a) 


or from the equation (B20b) 


P«($) - 


^ff ^ 1 

K eff^0 




1 


(C4) Kinetic distortion factors for the various area-ratio methods can be obtained 
from the computed quantities given in steps (Cl) and (C2). The procedure is as follows: 
(a) Determine NQ(r”) by using equation (B50b): 


N®(r) = f dE — <Pq(F, E) 
° J v(E) u 

(b) Determine N^(F) by using equation (B52a): 

Ng(r)= /dE-i-<pP(f,E) 
U J v(E) u 


(c) Determine the kinetic distortion factor for the Gozani and Garelis-Russell 
area-ratio methods based on fundamental mode values, by using equation (B53): 


,G0 


(r) = 


N^(r) 


00 


NS(r) F 


00 


z£V) = z^V) 

(d) Determine the kinetic distortion factor for the Sjostrand area-ratio method 
based on fundamental mode values, by using equation (B72’): 


GO/ 



(T) 


Ng(r) 

Ng(r) 
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The experimental procedure for obtaining reactivities from observed fundamental 
mode decay constants and measured prompt- and delayed-neutron fluxes is now de- 
scribed. The steps in the experimental procedure (El to E 6 ) are as follows: 

(El) Using a repetitively pulsed -neutron source and a v(E ) _1 neutron detector, ob- 
tain the detector response as a function of time. By means of a data processing code 
(ref. 13), the following information is obtained: 

(a) Fundamental prompt-mode decay constant, (agj 

(b) Prompt-neutron density for all modes, Np(r", t) 

(c) Delayed-neutron density for all modes, N^r", t) 

(d) Fundamental prompt-mode neutron density, Ng(F, t) 

(e) Parameter QqN^(F) defined by equations (B34), (B52a), and (B31a) and ob- 

tained experimentally by fitting the equation 


Ng(F,t) = Q q exp 



(B74) 


(E2) Compute an "experimental" value of the reactivity in dollars by using (oigj 

and computed values of /3 Q and A Q obtained as in step (Cl). Then, the "experimental" 
value of reactivity based on the inhour method is given by 



($) 



(B75) 


(E3) Compute the reactivity by using the Gozani extrapolated-area-ratio method. 
This reactivity is given by evaluating equation (B51), that is, 



($) 



Ng(F, t)dt 
N d (r,t)dt 


Z 



f“o) 

= - 7=5 [z?VS] (B76) 

Since the delayed-neutron distribution N d (F,t) is assumed to be constant with time, the 
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integral in the denominator of equation (B76) can be evaluated for a pulse repetition rate 
R. Then, this gives the reactivity as 



($) 


Q 0 N §(?) 



N d (r,t) 


z£°(F) 


R 


(B76') 


GO -*• 

The kinetic distortion factor Zq (r ) is defined in terms of computed quantities and is 
given in step (C4). 

(E4) Compute the reactivity by using the Garelis-Russell technique. This means 
finding a value of x such that the following equality is satisfied: 



N p(F,t)[exp(xt) - l]dt = N d (r,t)dt 


(B54) 


Since the delayed-neutron distribution N d (F, t) is assumed to be constant with time, the 
integral on the right side of equation (B54) can be evaluated for a pulse repetition rate R. 
This gives 


/ CO 

N p (r,t)[exp(xt) - l]dt = 


N d (r,t) 

R 


(B54') 


Then the reactivity in dollars is given by 

Z£ R (F) (B77) 

where Z^ R (F) = Z^°(r ). 

(E5) Compute the reactivity by using the Sjostrand area-ratio method (refs. 1 and 12). 
This reactivity is given as 


GR/a\ 
P s ($)= 


( a P\ 

{ a o) E 


+ 1 



($) =- 



• Z S0 (r) 


(B78) 


SO — ► 

where Z (r) is given by equation (B72). If N d (r,t) is constant with time, equation (B78) 
can be evaluated for a pulse repetition rate R. This gives 
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(B78*) 


P s ($) - 


-j N P (r - twt 
N d (r,t) 


Z S0 (F) 


R 


(E6) Compute the reactivity by using the method of Simmons and King (ref. 5). This 
reactivity is given by 


SK/ <t \ 

P s ($) = - 


[op\ 

™E 

(A 


+ 1 


L C 

E 


(B79) 


,DC 


where is the measured fundamental prompt-mode decay constant of the system 


'E 

at delayed critical. 
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